Vortex avalanches and magnetic flux fragmentation in superconductors 
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We report results of numerical simulations of non isothermal dendritic flux penetration in type-II 
superconductors. We propose a generic mechanism of dynamic branching of a propagating hotspot of 
a flux flow/normal state triggered by a local heat pulse. The branching occurs when the flux hotspot 
reflects from inhomogeneities or the boundary on which magnetization currents either vanish, or 
change direction. Then the hotspot undergoes a cascade of successive splittings, giving rise to a 
dissipative dendritic-type flux structure. This dynamic state eventually cools down, turning into a 
frozen multi-filamentary pattern of magnetization currents. 



PACS numbers: 74.20. De, 74.20. Hi, 74.60. -w 

The formation of a macroscopic current-carrying crit- 
ical state in type II superconductors occurs via pen- 
etration of the magnetic flux front of pinned vortices 
from the surface of the sample. This process is con- 
trolled by a highly nonlinear electric held-current den- 
sity (E — J) characteristics E(J,T, B), which together 
with the Maxwell equations determine macroscopic elec- 
trodynamics of superconductors. [H-SI Propagating mag- 
netic flux causes Joule heating, giving rise to global flux 
jumps and thermal quench instabilities which are crucial 
for stable operation of current-carrying superconductors 
Recent advances in the magneto-optical imaging 
have revealed a new class of instabilities of the critical 
state, including magnetic macroturbulence p,7|, kinetic 
front roughening 0], magnetic avalanches and dy- 
namic dendritic structures |K3-q3|. The latter have been 
observed both on high-T c (Y \Ba2Cu3O7 |lC|]) and low- 
T c (Nb |ll|), superconducting films, and most recently 
on the newly discovered MgB2 p2| . These instabili- 
ties are rather characteristic of superconductors, besides 
they also display remarkable similarities with other den- 
dritic instabilities in crystal growt h [Pi , nonequilibrium 
chemical and biological systems [Q, and crack propa- 
gation ifnU . This analogy is not surprising, because the 
non isothermal flux dynamics in superconductors is de- 
scribed by two nonlinear equations for the magnetic in- 
duction B(r,t) and temperature T(r,t) somewhat simi- 
lar to those describing dendritic structures in reaction- 
diffusion systems |Q. 

We performed numerical simulations of coupled equa- 
tions for the magnetic induction B(r, t) and temperature 
T(r,t) and found a new mechanism of flux fragmenta- 
tion in superconductors, which is different from the well- 
known bending instability of moving interface between 
two phases [ fl3[ , but results from the generic distribution 
of magnetization currents in the critical state. Our re- 
sults give insight into unstable flux penetration, vortex 
micro avalanches and flux jumps in superconductors. 

The effect is best illustrated by considering a slab in 
a magnetic field Bq (Fig. 1), for which distributions of 



the z-component of magnetic induction, B(r, t), and tem- 
perature T(r, t) are described by the Maxwell equation 
coupled to the heat diffusion equation. 

Cd t T = divKVT - h(T - T )/d + JE{ J, T), (1) 
d t B = -cV x E(J, T), J = (c/4tt)z x VS. (2) 

Here C(T) is the heat capacity, k(T) is the thermal con- 
ductivity, h{T) is the heat transfer coefficient to the 
coolant held at the temperature T , d — A/P, A is 
the area of the sample cross section, P is the perime- 
ter of the cooled sample surface, B is the z-component 
of the magnetic induction, and E(J,T) is the modulus 
of the electric field, which essentially depends on both 
the local temperature T(r, t) and the current density, 
J(r,t) = (J 2 + Jy) 1 / 2 - The spatial derivatives are taken 
with respect to x and y, while the term h(T — To)/d ac- 
counts for the surface cooling Jig] , We will consider the 
case of high magnetic fields Bo, much greater then the 
field of full magnetic flux penetration, B p ~ A-nJ c d/c, for 
which E(J, T, B) re E(J, T, B ). 

The evolution of T(r,t) and B(r,t) is mostly deter- 
mined by the E(T, J) characteristic, which accounts for 
high-resistive flux flow state at J > J c and low-resistive 
flux creep state at J < J c , where J c is the critical cur- 
rent density. The effects considered in this work are not 
very sensitive to the details of E(J,T), so for the nu- 
merical simulations of Eqs. (Q)-^) we take the following 
interpolation formula expressed in terms of observable 
parameters J c , J±, and p: 



E = pJi ln[l + cxp(J - J c )/Ji] 



(3) 



Here Ji(T) = dJ/dlnT is the dynamic flux creep rate, 
and p(T) = p n B/B C 2 is the flux flow resistivity. Below 
the irreversibility field B < B* , where J\ ^ J c , Eq. 
(^) reproduces the main features of E( J, T) observed in 
experiment, giving a linear flux flow dependence E = 
(J — J c )p for J > J c and the exponential dependence 
E = E c exp(J - J c )/Ji for J < J c . 

The similarity of Eqs. (|)-(|) with generic reaction- 
diffusion equations |Q, is due to thermal bistability of 
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superconductors g, for which the heat balance condi- 
tion (T - T )h/d = JE(T,J) in the right-hand side of 
Eq. (0) is satisfied for 3 different temperatures T, as 
shown in Fig. f. Here the points and 3 correspond 
to two stable uniforms states: a cold superconducting 
state with T ks To and a hot flux flow/normal state T3 
self-sustained by Joule heating. As seen from Fig. 1, 
the bistability occurs, if the current density J exceeds 
a threshold value J m , for which pj^ ~ h(T* — To)/d. 
Hence, J m ~ [(T* - T )h/ 'dp} 1 / 2 , where T*(B) is the 
irreversibility temperature at which J C {T*) = 0. The su- 
perconducting state is unstable with respect to the hot 
spot formation, if a s — (J c /J m ) 2 > 1. For typical pa- 
rameters of HTS films at T = 4.2K: p ~ 100/ificm, J c = 
10 6 - 10 7 A/cm 2 , h ~ lW/cm 2 K, T* - T Q ~ 50 - 100K 
and d = lp,m, we obtain a s ~ 10 2 — 10 4 , thus the ther- 
mal bistability is a characteristic feature of both HTS and 
LTS, especially films because of their higher J c values ||. 

We consider the case of weak Joule heating, T(x, y, t) ~ 
T c , for which we take into account only the most essential 
temperature dependence of E(T), while C(T), h(T), and 
k(T) can be taken at T = T . Then Eqs. @)-<g) can be 
written in the following dimensionless form 



Tb = d x [r(j,e)d x b] 

e = v 2 e- 



1 + aj 2 r{j,9). 



(4) 
(5) 



Here 6 

\{d x b) 2 



= (T - T )/(T* - T ), b = B/B t and j = 
- {dyb) 2 ] 1 / 2 are the dimensionless temperature, 
magnetic field, and current density, respectively, and 
B t = 4irJiL h /c. The derivatives in Eqs. (g) and @ 
are taken with respect to normalized time t/th and co- 
ordinates x/Lh and y/Lh measured in the thermal units 
t h = Cd/h and L h = (dn/h) 1 / 2 . The evolution of 0(r,t) 
and fr(r, t) is controlled by two dimensionless parameters: 

_ 4ttk pJ\d 
T ~ P Cc 2 ' ° ~ h{T* - T ) ( ' 

Here r is the ratio of the diffusivities of heat and mag- 
netic flux, and a quantifies Joule dissipation. For Nb 
films of Ref. jnj, we obtain r ~ 10 at 4.2K, with r 
rapidly decreasing with increasing T . For HTS at 77K, 
we get t ~ 10~ 4 — 10~ 5 <C 1. The nonlinear resistivity 
r(j, 9) = ln[l + exp(j — j c (&))]/j obtained from Eq. (3) 
has asymptotics r = 1 —j c /j in the flux flow (J > j c ) and 
r = cxp(j — j c )/j in the flux creep (j < j c ) states, where 
jc = Jc(T)/J\. We linearize j c = jo(l — 9) around To, 
neglecting the temperature dependencies of J\ , p jl7j . 

We performed 2D numerical simulations of Eqs. (g) 
and (||) to calculate propagation of a magnetic hotspot 
of resistive phase across a superconductor (see Fig. 1). 
The process is initiated by a local heat pulse applied to 
the sample surface, which models the experiment by Lei- 



derer et al., |10|, who triggered the magnetic dendrite 
instability by a laser pulse. To address the effect of ma- 
terial inhomogeneities, we considered both uniform su- 
perconductors with J c independent of spatial coordinates 



and nonuniform superconductors with the critical current 
density periodically modulated over macroscopic scales 
2ir/k, much larger than the spacing between flux lines, 
J c (x,y,T) = J c0 {T)[l + esm{kx)sm{ky)} withe < 1. The 
latter case also models superconducting films with peri- 
odic arrays of holes, which have recently attracted much 
interest Jl^,[ll| . The resulting evolutions of the tempera- 
ture distributions shown in Figs. 2 and 3 display a rather 
striking behavior described below. 

For t <C 1, the heat pulse applied to a uniform su- 
perconductor triggers a hot spot propagation across the 
sample, as shown in Fig. 2a. Such propagation ap- 
pears stable until the hot domain reaches the center of 
the sample, where magnetization currents change direc- 
tion. Then the resistive domain undergoes a cascade of 
successive splittings into alternating stripes of low and 
high electric fields and temperatures. After each succes- 
sive splitting, the part of the resistive domain near the 
central line cools down, but then the hot filaments of re- 
sistive state start propagating again from the upper part 
of the resistive domain through the preceding dendritic 
structure toward the central line. However, each time the 
hot filaments cross the central line, they split again, caus- 
ing new dendrites of alternating low and high-J filaments 
to grow as shown in Figures 2b and 2c. Eventually the 
Joule dissipation causes the electric field in the hot den- 
dritic structures to decay below the threshold, at which 
further hot spot propagation becomes impossible, thus a 
frozen entangled pattern of current filaments forms. 

As follows from Figures 2d-2f, spatial inhomogeneities 
in J c can bring about new features of the hot spot propa- 
gation, causing fragmentation of the hot spot even before 
it crosses the central line. The periodic modulation of 
J c (x,y) also gives rise to additional side branching and 
preferential flux propagation along the nearest neighbor 
directions at angles ±45°, causing further interconnec- 
tion of neighboring hot spot branches. This effect is sim- 
ilar to that observed by Vlasko-Vlasov et al on supercon- 
ducting films with periodic arrays of holes and magnetic 
dots ]l9|[ . For r <C 1, the diffusion of magnetic flux in 
the normal state occurs much faster than the heat dif- 
fusion. Therefore, magnetic flux rapidly penetrates the 
hot regions of the filamentary current structure in Fig. 2, 
forming dendritic flux front patterns reminiscent of those 
observed in magneto-optical experiments. ||Io|-|l2|| 

The flux fragmentation can be described as follows. 
For r <C 1, the electric field E becomes nearly poten- 
tial, V x E rj 0, thus the magnetization currents tend 
to bypass propagating hot resistive domain, similar to 
that of a steady-state current flow past a planar obstacle 
p0| |. Thus, the current density in the resistive domain 
decreases, forcing the excess magnetization currents to 
flow along the domain interface. The high interface cur- 
rents cause strong local enhancement of electric field and 
dissipation, widening the resistive domain near its end 
accelerating the propagation velocity. At the same time, 
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the temperature in the center of the resistive domain de- 
creases, facilitating recovery of the superconducting state 
when the domain crosses the central line, where magne- 
tization currents change direction. In this region the in- 
terface currents at the bottom part of the domain are 
partly compensated by the opposite magnetization cur- 
rents, which strongly reduces Joule dissipation and stops 
hot spot propagation. As a result, a triangular region at 
the tip of the resistive domain becomes superconducting, 
and then the process repeats as described above. 

For t ^> 1, the hot spot propagation occurs at the 
frozen distribution of magnetic fields and currents, and 
the hot spot dynamics changes as shown in Fig. 3. In- 
stead of the multiple splitting characteristic of r <C 1, 
the heat pulse at r > 1 first initiates stable propagation 
of a resistive hot spot, then its splitting into two parts 
which move apart and eventually disappear. This behav- 
ior occurs if the energy of the heat pulse Q is below the 
critical value Q c . For Q > Q c , the heat pulse creates a 
larger hotspot which then expands and propagates over 
the entire sample ||. 

The dendritic flux penetration described above can 
also be regarded as a micro avalanche of a macroscopic 
bundle of vortices, which does not trigger a global flux 
jump instability or thermal quench of the whole sample. 
Such micro avalanches only cause local transient tem- 
perature spikes, leaving behind the frozen flux dendrite 
structures shown in Figure 4. Each micro avalanche thus 
results in a partial flux penetration, which reduces the 
total magnetic moment of a sample and manifests itself 
in steps on magnetization curves M(B). Notice that mi- 
cro avalanches in increasing magnetic field may also be 
initiated by surface defects (regions with lower J c ), which 
are sources of excess steady-state Joule dissipation. Such 
common defects, which have been revealed by magneto- 
optical imaging of HTS Q , can trigger both global flux 
jumps and local vortex micro avalanches. The micro 
avalanches cause steps in M(B) as observed on Nb films 
at low temperatures Q . After many microavalanches the 
critical state eventually turns into a frozen "turbulent" 
current structure built of individual dendritic fragments, 
like that in Fig. 4 

The results this work may capture the essential physics 
of dendritic flux instability observed on Y BazCu-sOi jlO| , 
Nb [JTTJ] and MgB 2 @ films, although, for a more quan- 
titative comparison of the theory with experiment, other 
factors should also be taken into account. Indeed, the 
experiments [^0|-^2[ correspond to thin films in low per- 
pendicular magnetic fields, which requires the account of 
the nonlocal flux diffusion ||] and the geometrical barrier 
p2[ , whereas our model describes a slab in high paral- 
lel magnetic fields unaccessible to magneto-optical tech- 
nique. Another intriguing result discovered in Ref. JTcj ] 
is a superfast dendritic flux propagation with the veloc- 
ities exceeding the speed of sound c s . This also requires 
invoking additional mechanisms, because the thermal ve- 



locities of hotspot propagation!; ~ J[np/(T* — Tq)] 1 / 2 /C 
H is smaller than c s . The superfast dendritic flux prop- 
agation might be due to electron overheating, so that an 
equation similar to Eq. (1) actually describes the elec- 
tron temperature, higher than the lattice temperature 
To. Then the lattice heat capacity C in the above esti- 
mate for v( J) is to be replaced by a much smaller electron 
heat capacity, which increases v(J) by 1-2 orders of mag- 
nitude. This situation may occur at low T, if the time 
of the electron-phonon energy relaxation becomes larger 
than the thermal time th- p3fl . 

In conclusion, we proposed a new mechanism of mag- 
netic flux fragmentation in superconductors. The insta- 
bility manifests itself as a vortex micro avalanche, accom- 
panied by a transient local Joule dissipation and eventu- 
ally resulting in a frozen multi-filamentary structure of 
magnetization currents. These effects give rise to den- 
dritic flux penetration into superconductors, partial flux 
jumps and steps on magnetization curves. 
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FIG. 1. Graphic solution of the heat balance condition 
(T — To)h/d — JE(6,J) in a superconductor. The power 
of the Joule heat release JE(J, 9) is plotted as a function of T 
for Ji < J m and J2 > J m , where E(J,8) is given by Eq. (|3|), 
and 6 = (T — Tb)/(T* —To). Inset shows the sample geome- 
try, magnetic field is parallel to z axis, a;-axis is directed along 
the sample surface and y is perpendicular to the sample sur- 
face. The local heat pulse (indicated as big arrow) triggers the 
magnetic hotspot propagation across magnetization currents 
presented in detail in Figs. 2 and 3. 
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2. Gray-coded dynamic temperature maps 9(r 



the flux fragmentation instability in homogeneous system for 
t/th =8 (a), 20 (b), 34 (c) at a = 0.008, r _1 = 150 and 
the initial current in the Bean state, jo = 18 (white corre- 
sponds to > 1, and black to 9 — 0). Domain of integration: 
— 1501/h < y < 1501/h, < x < GOOLh, periodic boundary 
conditions in x direction; no-flux d9/dy at y = ±150Lft, for 
the temperature and b = const for the magnetic field. Each 
panel shows the upper half of the sample and 1/3 of total 
lenght: < y < 150L h , 200L h < x < 400L h . The magneti- 
zation currents change direction on the bottom part of each 
panel, i.e for y = 0. Same for the system with periodic mod- 
ulation in J c (x,y) for 27r/fc = 30L h , e — 0.5, t =6(d), 16(e), 
28(f) and j = 20. 
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FIG. 3. Same dynamic temperature maps as in Figure 2 
for t = 1, a = 0.006, j = 18, and t/t h =8 (a), 20 (b), 34 (c). 
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FIG. 4. A typical "frozen" current configuration after com- 
pletion of the fragmentation for r _1 = 600. The gray shades 
change from black (J = 0) to white (maximum J). Full do- 
main of integration is shown. Black line in the middle is due 
to reverse of direction of the magnetization current. 



